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Abstract. In contrast to the field, the binaries in dense stellar systems are frequently not primordial, and could be either 
dynamically formed or significantly altered from their primordial states. Destruction and formation of binaries occur in 
parallel all the time. The destruction, which constantly removes soft binaries from a binary pool, works as an energy sink 
and could be a reason for cluster entering the binary-burning phase. The true binary fraction is greater than observed, as a 
result, the observable binary fraction evolves differently from the predictions. Combined measurements of binary fractions 
in globular clusters suggest that most of the clusters are still core-contracting. The formation, on other hand, affects most 
the more evolutionary advanced stars, which significantly enhances the population of X-ray sources in globular clusters. The 
formation of binaries with a compact objects proceeds mainly through physical collisions, binary-binary and single-binary 
encounters; however, it is the dynamical formation of triples and multiple encounters that principally determine whether the 
formed binary will become an X-ray source. 
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BINARIES EVOLUTION: DESTRUCTION 

For low density environments, like solar neighborhood and for solar-type stars, the binary fraction is known to be 
^ 50% (e.g., [1]). Similarly, open clusters show relatively large measured binary fractions; as large as 70% in some 
open clusters; a clear anti-correlation with the age of the open cluster has been found [2]. Most globular clusters have 
typically significantly lower measured binary fraction than open clusters, although for some it is comparable to that of 
open clusters - e.g., in a core-collapsed NGC 6752, the upper limit for the binary fraction can be as high as 38% [3] 
(although see also [4]); in a young and sparse globular cluster in Ter 7, the binary fraction is estimated to be 51% [5]. 
For globular clusters, the measured binary fraction was found to anti-correlate with the total cluster luminosity [6]; 
the correlation with the cluster collisional parameter was found to be only marginal. All together, measurements of 
binary fractions in open and globular clusters suggest that binaries in clusters deplete with time and it happens more 
efficiently in more massive clusters. 

There is no theoretical study that would suggest that the initial binary fraction should be significantly different at 
the moment of star formation between proto-open and proto-globular clusters ' . The consistency between the initially 
large binary fraction and currently measured low binary fraction can be explained by the efficient binary destruction 
that continues throughout the entire cluster evolution. 

Indeed, let us consider a population of primordial binaries with the initial periods P = 0. 1 — 10^ and flat mass ratio 
distribution, evolved for 10 Gyr using StarTrack [8]. By this age, significant fraction of initial binaries is destroyed 
through evolution, mainly due to mergers, or, in case of massive stars, in supernovae (see Fig. 1). If we define the 
hardness of a binary 77 as the ratio of the binary binding energy to the kinetic energy that a O.5M0 object moving 
with V = \/3 X 10 km/s has (this energy is about an average kinetic energy of an object in a typical dense cluster 
and is referred here to as kT), then it can be seen that by the age of 10 Gyr, almost all the hardest (rj > \QQkT) 
binaries are destroyed via evolution. Now we place a population of these binaries - with primary masses drawn in 
accordance to initial mass function (IMF) from [9], and with their initial eccentricities distributed thermally - in a 
confined volume and keep the concentration of objects constant «c = 10^ pc^^, constant velocity dispersion of 10 
km/s and disallow the objects to escape form this test volume. If we let the same stars evolve for the same 10 Gyr, but 
allow them to interact dynamically, we find that - not surprisingly - almost all the soft binaries (tj > IkT) are destroyed 



' Recent studies showed that there could be more than one period of star formation in globular clusters, and different generations could have been 
fornied with different binary fractions [7]. 
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FIGURE 1. Probability density for a primordial binary to survive 10 Gyr as a binary: the effect of the evolution only (left), the 
effect of dynamics (right). Solid lines shows rj, binary hardness, and dashed lines show Tcoih collision time. 

by dynamical encounters, except those with the colhsion time^ Tcou of the same order as the cluster age (for the Fig. 1, 
for Tcoii, an average object is assumed to have mass of O.SM©). Roughly, the binaries with Tcou more than few Gyr, 
are not perturbed by dynamical encounters, and with Tcou less than few Gyr - only the hard binaries - did undergo a 
dynamical encounter but survived. 

Clearly the efficiency of the destruction of binaries in reaUstic dense stellar systems is somewhere in between these 
limiting examples, and also strongly affected by the evaporation of both binaries and single stars from the core due to 
high post-encounter velocities or supernova kicks, and by continuing segregation of binaries and single stars from less 
dense halo, as will be discussed below. 



Evolution of hard binaries 

Two main approaches currently used to model the dynamics of a dense stellar system are direct A^-body, the most 
widely known realization is the family of NBODY codes written by S. Aarseth^, with the most recent version A''B0Z)F6; 
and Monte Carlo (MC) methods, where encounters itself are usually treated with a direct few-body integrator (e.g., 
[10, 11]). MC methods are more capable of treating large populations and have been already used to make models of 
observed globular clusters (e.g., M4 [12], NGC 6397 [13], 47 Tuc 1 14]), while direct N-body is usually used to model 
open clusters (e.g. NGC 188 [15], M76 [16]), but is starting to be also apphed to globular cluster (e.g. NGC 6254 
[17]). Here, we will refer MC to the specific realizations that self-consistently model the global evolution of a cluster, 
using MC techniques to sample the stellar distribution function when applying the effects of two-body relaxation [18]. 
This not to be confused with a simplified Monte Carlo method with a two-zone model of a cluster as in [19]. In all 
the cases, the dynamical codes assume as well the use of a (binary) stellar evolution (population synthesis) code of a 
different complexity. 

Due to still persisting computational power limitation, even with nowadays super-computers, including those 
equipped with GPUs, the numerical A^-body models of globular clusters are traditionally calculated either using 
small number of primordial binaries or small mass clusters (up to 100,000 stars, e.g. [20]), or rather very sparse. The 
limitation so far comes not from the stellar evolution, which is either parametrized or tabulated, but from dynamics. 



^ this is the time-scale for a binary to undergo a strong encounter with another object 
' See http://www.ast.cam.ac.uk/'^sverre/web/pages/nbody.htm 
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The maximum initial binary fraction used for a direct A^-body modeling is 50% for an (open) cluster made of 12,000 
single and 12,000 hard binary stars and with core density nc = 100 — 350 pc^-' [21]. For relatively more massive 
clusters, up to 100,000 stars (a factor of 10 short compared to observed massive and dense clusters) initial binary 
fraction is usually 5%. (e.g. [20]). 

With the same initial conditions for hard binaries (low or medium density clusters, binary fraction < 50% ), MC and 
direct A^-body do agree each other qualitatively and quantitatively: the core binary fraction increases with time (see 
the comparison made in [18]). For larger initial binary fractions, MC shows the same behavior for almost all initial 
binary fractions up to = 90% (here, all primordial binaries are hard). The simulations were performed for cluster 
models with initially 10^ objects and with IMF extends from 0.15 to only 18.5 Mq (MC is not designed to handle 
subsystems of rare objects, like black holes (BHs), which could be produce by stars more massive than I8.5M0); the 
highest initial core density lO'*^ pc^^. 

Interesting enough to note that the increase of the core binary fraction with time was due to the two reasons. The 
first was theoretically expected to play the most significant role - the segregation of binaries from the halo to the core. 
However, especially in the case with binary fraction no being close to one, the main net effect for the binary increase 
was at large due single stars evaporation from the core (in fact, just a few % of original core single stars remain in the 
core, the rest is evaporate either to the the clusters outskirts, or lost from the cluster completely, then for binaries about 
a half of original hard binaries remains in the core). By the end of the core mass was decreased by at least an order of 
magnitude with time. 



Comparison with observations: is there a steady binary-burning phase? 

In cases when MC models entered the binary-burning phase (not a common results among all models) before a 
Hubble time, it was found that the core binary fraction in this phase steadily decreases with time. This behavior is 
consistent with the results of simplified two-zone Monte Carlo model described in [19], where the binary-burning 
phase was imposed via adopted constant with time cluster core properties. It is hkely then that efficient the hard binary 
fraction depletion can proceed only in a post-collapsed cluster. How conmion can this be? 

As it was mentioned, in most of numerical simulations with hard binaries, the steady binary-burning regime is not 
observed. A cluster evolution towards the core collapse can take up to a Hubble time. One other hand, the theoretical 
predictions for a cluster in a binary-burning phase do not match well with observations. On Fig. 2 we show the 
comparison of the cluster cores predicted by theory for clusters in a binary-burning mode and of the observed cluster 
cores. It can be seen, that the theoretical cluster cores are significantly (by an order of magnitude) are smaller than the 
observed ones. This likely indicates that the most of globular clusters in Milky Way is not in a binary-burning phase 
[22]. 



Soft binaries: can they still be neglected? 

The described above results were obtained using only initially hard binaries, as soft binaries have been traditionally 
ignored in both Monte Carlo and direct A^-body simulations. Mainly, this is because soft binaries are expected to be 
destroyed on a very short time-scale and as such are not expected play a role on the globular cluster dynamics as a 
whole. The other reason, specific for A^-body, is that each wide (soft) binary requires enormous calculation power, 
significantly slowing down the cluster's simulation. 

However, largely unexpected result has been obtained in the model that included soft binaries, and had an initial 
binary fraction (pg = 90%. In this case, after the rapid initial contraction of the core, the cluster fairly quickly (compared 
to the models with only hard initial binaries, on the scale of just a few Myr) enters into a long-lived binary-burning 
phase keeping almost constant binary fraction of 40% [18]. The results were obtained for a cluster model with 5x10^ 
stars. 

The physical reason for a different behavior of a cluster with a substantial contain of initial soft binaries is that 
this soft component acts, while present, as a significant cluster energy sink. Soft binaries indeed are ionized quickly 
through encounters, absorbing clusters' energy. As such, if the core was bom with a significant number of soft binaries, 
it will rapidly contract on the time-scale of soft binaries destruction. 
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FIGURE 2. The comparison of the theoretical relation for the cluster cores at binary-burning phase and the observed cluster cores; 
(j>S is the core binary fraction. Credit for this figure is to John Fregeau, with thank to Craig Heinke for compiling the observational 
data with the same assumption (e.g., flat mass ratio) ([4, 19, 23, 24, 25, 26, 27, 28, 29]). The theory curve comes from [30]. 

What is a true binary fraction? 

It is not widely recognized, but theoretically obtained binary fractions are very different from the observed ones. 
When describing simulations and its match to observed globular clusters, therefore, it is important to distinguish 
between a "true" theoretical binary fraction - how many binary objects are present in a simulation, independently on 
their characteristics, such as mass and luminosity, and an "observed" binary fraction - how many binaries would be 
found by an observer if this simulation would be a real cluster. 

It can go either directions. For instance, when observers find binaries in a specific period range, an overestimation 
of binaries fraction can be done if it is assumed that the birth period distribution has been conserved for hard binaries, 
and only soft primordial binaries are destroyed. In [19] it was shown that even hard binaries get destroyed, when wider 
hard binaries are destroyed by dynamical encounters (a significant fraction of binaries up to IQkT is very vulnerable), 
and the evolution takes away very hard or intermediately hard but massive binaries (see also Fig.l). E.g., in [23], the 
observationally derived binary fraction for 47 Tuc was found to be 13%, but if the change in periods distribution is 
taken into account, then the extrapolation from the observed binaries can give only the binary fraction of about 6%. 

On the other hand, [18] had demonstrated that even in simulations with only hard binaries, when the overall "true" 
binary fraction in the core increases, the "observed" binary fraction decreases for all simulations with the initial hard 
binary fraction above 40%. For instance, a medium density cluster with initial core density 10^^ pc^^ and initial binary 
fraction 90%, evolved to only 30% of observable binaries in the core before the clusters was tidally disrupted. The 
"true" binary fraction though in this case approached almost 100%. "Hiding" binaries here include binaries where two 
main sequence (MS) stars with a low-mass ratio would blend for an observer into one MS star, or binaries with one or 
two dim compact object as a companion. 

BINARIES EVOLUTION: FORMATION OF COMPACT BINARIES 

Although most dynamical interactions in dense cluster cores tend to destroy binaries, some can form binaries from 
singles stars, and many modify them. The most important processes are those that lead to a formation of binaries both 
have a compact companion and are compact enough to start mass transfer. Such binaries can be detected not only in 
our MUky Way, but in globular clusters in distant galaxies as well, giving us a link between the statistical properties 
of internal parameter of globular clusters and the efficiency of dynamical encounters. In particular, low-mass X-ray 
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binaries (LMXBs) are formed in globular clusters at the rate of 100 times exceeding that of in the field, per stellar mass 
unit [31]. The most of binary millisecond pulsars (bMSPs) - the likely termini of LMXBs evolution - detected so far 
are located in globular clusters [32, 33, 34]. It was even proposed that all LMXBs and bMSPs are formed in globular 
clusters [35]. The importance of dynamical encounters for LMXBs formation seems to be well estabUshed, as the 
number of LMXBs was found to correlate well with the cluster dynamical properties [36] for all non-core-collapsed 
clusters. [37] explained the exception of core-collapsed clusters by our misunderstanding of clusters current dynamical 
states and suggested that most globular clusters are still in the process of core contraction. More recent observations 
showed that the overabundance of LMXBs in core-collapsed clusters is statistically significant, and the number of 
X-ray sources in such clusters is almost independent on the cluster's dynamical properties [38]. The only LMXB with 
a BH and a white dwarf (WD) was detected in a globular cluster [39], although their theoretically predicted numbers 
in the field should vastly exceed the number of LMXBs with a MS companion [40, 41], and 17 LMXBs with a MS 
companion in the Milky Way are observed. 

There are several processes that lead to close binary formation: (1) binary companion exchange; (2) physical 
colhsion; (3) tidal capture and (4) three-body binary formation. The evolution of a dynamically formed binary can 
be further perturbed by dynamical encounters: a binary can be hardened and its eccentricity can be "pumped" via 
multiple non-strong encounters; it also can become a member of a triple. The top channel to form a binary with a 
neutron star (NS) or a WD is a binary companion exchange, if binary fraction is at least few % [42, 43]. A tidal 
capture (TC) operates only within a narrow range of periastra during a single-single star encounter [42]; as such, 
this channel does not account for more than a few % of all formed binaries with a compact companion. It has been 
argued that TCs could play an important role in the formation of ultra-luminous X-ray sources, with an intermediate- 
mass BH (IMBH) acquiring a companion through a TC [44, 45]. However, in most of TCs events with an IMBH, 
the energy dissipation rate in a captured star greatly exceeds the star's Eddington luminosity [45]. The fate of such a 
binary depends then on where inside the star the energy was deposited, and could be a merger [46]. Three-body binary 
formation is capable of creating a hard binary, but the formation rate of binaries which are hard enough to become 
X-ray sources is small, even for massive black holes [47]. Hardening, through multiple encounters, can shrink a hard 
binary towards it Roche lobe overflow, though a fraction of the binaries that can successfully survive this path is not 
very large [47]. In the following we will review in more detail the processes that received attention in recent years: 
formation via physical colhsions and effect of the triples formation on compact binaries. 



Formation of close binaries via physical collisions 

The formation of binaries via a physical collision between a red giant and a compact object has been proposed first 
by [48], in order to explain ultra-compact X-ray binaries (UCXBs) formation. In this scenario, a colhsion leads to the 
formation of a bound system that later might experience a common envelope and form a tight binary (note though that 
for low-mass giant remaining in globular clusters at current age, mass ratio favors a dynamically stable mass transfer). 
Smooth particle hydrodynamics (SPH) simulations showed later that the stellar envelope can be disrupted or fully 
removed in close encounters, with an eccentric binary formed as a result [49, 50]. Later, considering subgiants and 
early giants obtained using a stellar evolution code, it was shown that (SPH) physical collisions between them and 
NSs, with small impact parameter rp < 1.5/?rg where Rrq is a giant radius, always lead to the complete expulsion of 
the envelope during a colhsion, resulting in a tight and eccentric NS-WD binary; this binary shortly thereafter decays 
its orbit and starts mass transfer [51, 52]. Formation rate of UCXBs by these encounters is sufficient to explain the 
observed number of UCXBs in Galactic globular clusters or LMXBs in globular clusters near other galaxies [51]. The 
formation rates for UCXBs, being consistent with the observed number, at the same time predict the number of bMSPs 
well exceeding the observed number of bMSPs in globular clusters, and it was suggested that this type of LMXBs does 
not produce a radio bMSP [43]. 

When the first BH-WD LMXBs was detected, it was plausible to consider a similar formation mechanism. Detailed 
SPH numerical simulation showed that even though in most of encounters with rp up to ~ 5/?rg a bound system 
is formed, only a small fraction of them (with rp up to ~ /?rg) will lead to a formation of binaries that are tight 
enough and eccentric enough to start mass transfer in isolation [47]. A formation rate of BH-WD LMXBs through 
only physical collisions is about order of magnitude lower then required to explain the observed formation rate; the 
binary exchange channel does not provide binaries compact enough to start the MT in isolation, so a sequence of 
post-coUisional dynamical encounters is necessary. 
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Role of triples in compact binaries formation 



One of the outcomes of a binary-binary encounter is the formation of a triple system (e.g. [53]), and some of the 
formed triples will be hierarchically-stable [54]. In a typical dense cluster with central density 10^ pc~^ and binary 
fraction of ^ 10%, at 10 Gyr, about 5% of all core binaries would have successfully formed a hierarchically stable triple 
during 1 Gyr [55]; over whole hfe of the cluster, about half of binaries could become a member of a triple. A typical 
formed triple has large ratio of orbital periods, Pout/Pin ~ 1000 and very large outer eccentricity, gout = 0.95 ± 0.05. 
About a half of formed triples is hard. 

The effect of triples formation on the evolution of binaries population in globular clusters has just started to be 
recognized over the last several years. Driving force behind the triples' effect on binaries is that, if a triple has a large 
enough inclination, its secular evolution is affected by Kozai mechanism [56]. This mechanism causes large variations 
in the eccentricity and inclination of the star orbits and could drive the inner binary of the triple system to Roche lobe 
overflow, when as result of the starting mass transfer a binary either merges before the next interaction, or start stable 
mass transfer. If one component of the inner binary is a non-degenerate star, tidal interactions can be important during 
the periastra, when Kozai-induced eccentricity is at its maximum. Shrinkage of a binary by a combination of Kozai 
cycling and tidal friction (KCTF; [57, 58, 59, 60]) can be responsible for production of short-period active stars, e.g. 
BY-Dra-type binaries [61]. KCTF would operate in a cluster if formed triple's Kozai time-scale [62] will be smaller 
than Tcou. As triples formation provides rather uniform space distribution for incUnations, about 1/3 of all the formed 
triples affected by Kozai mechanism. 

Blue stragglers, which are exotic stars with masses up to three cluster's turn-off mass but presumably still burning 
hydrogen, are detected both in globular and open clusters. They are beUeved to be either produced by collisions 
between stars in clusters or stable or unstable mass transfer between the components of primordial short-period 
binaries. Using a simplified Monte Carlo method, it was found that a fraction of blue stragglers which can be provided 
by KCTF is comparable to that from physical collisions and primordial binaries [55]. KCTF also has been proposed 
to have a dominant role in the formation of blue stragglers in open clusters and in the field [63]. Specifically, it 
could naturally explain the large binary fraction of long-period blue stragglers binaries; these results are based on an 
assumption of primordial population of triples rather than on the population of dynamically formed triples. 

As a star evolves, the role of the triple in the formation of a binary with such a star is increasing, and the probability 
that an inner binary would have a compact star as a companion is higher than for overall binary population [55]. 
E.g., if KCTF in a triple with an inner binary consisting of a WD and a MS star leads not to a merger but to a stable 
mass transfer, the number of formed and present cataclysmic variables can be increased by up to 50%, if compared 
to simulations where KCTF is not taken into account. If KCTF is taken into account, but leads to a merger, then the 
number of CVs could be depleted by third. 

In a 47 Tuc type cluster, 5% of all binaries with a NS companion would become a member of a hierarchically stable 
triple in a Gyr; in more dense cluster, Ter 5 - type, 15% of all binaries with a NS would be a member of a triple. In 
simulations, it appear that ^ 50% of all bMSPs were in triples at some point in their past [43]. 

In simulations of globular clusters, since there is as yet no "triple" population synthesis, a dynamically formed triple 
is simply broken, and the inner binary is shrunk to the minimum periastron it could have through simple Kozai cycle, 
if Kozai time scale is shorter than Tcou- So in current clusters simulations, the mass transfer could be induced by triple 
formation, but not forced by being continuously in a triple. 

There is, however, an UCXB 4U 1820-303, located in the globular cluster NGC 6624, with the orbital period P of 
685s [64, 65]. The stabihty value P/P = -(3.5 ± 1.5)10"^ yr"' makes certain that 685s is the orbital period [66]. 
Secondary star is a He WD 0.06 — 0.08 M© [67]. Note that standard scenario of Roche lobe mass transfer for a such 
otherwise rather usual UCXB binary implies positive period derivative. [66] found that, in addition, 4U 1820-303 has 
the luminosity variation by a factor of 2 at a super-orbital period Pso ~ 17 Id. It has been suggest that, since X-ray busts 
take place only at the flux minimum, the observed variability is due to intrinsic luminosity/accretion rate changes and 
not obscuration or changes of the projected area of the source due to precession. The ratio between super-orbital and 
orbital periods ( 22000) is too high for any kind of the disk precession at the mass ratio of the system [68, 69]. 

It has been suggested that this binary is part of a hierarchical triple [66], with the third body having a mass 
< O.5MM0 (based on the lack of its optical detection) and with the third body orbital period Pout ^ 1 . Id (as eccentricity 
modulation is expected to be on P^o — PonJP)- [70] studied this system in detail with the inclusion of the following: 
perturbation from a third body on a longer period orbit; the quadrupolar distortion of stars due to their intrinsic spins 
and the further quadrupolar distortion due to their mutual gravity; tidal friction in the equilibrium tide approximation; 
general relativity; mass transfer and gravitational radiation. They found that the long period naturally arises if the 
system is Ubrating around the stable fixed point in a Kozai resonance. The observed system including the long period 
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for eccentricity modulation can be reproduced with a NS with mass of 1 AMq, WD mass of O.O67M0 and third body 
mass of O.55M0. The semi-major axes of the inner binary is a = 1 .32 x lO'^cm, of the outer binary flout = 6.52a and 
corresponding period of /out = 0.1 Id. Initial eccentricity eo = 0.002, initial outer binary eccentricity Cout.o = 0.0001 
and initial inclination 65.66°. The amplitude of the eccentricity oscillations was found to be of ~ 3 x 10~^, which is 
sufficient to enhance mass transfer enough to produce the observed luminosity oscillations by factor of 2 [71]. 

Although there is not yet a detected triple with a BH as a companion of the inner binary, the role of dynamical triple 
formation for compact binaries with a BH could be even higher than for compact binaries with a NS companion. For 
binaries with similar masses and orbital separation ai and 02, the fraction of binary-binary encounters that result in a 
hierarchically stable triple formation is about 25% for all encounters that have a periastron within ^ 7(fli +02) (for 
details see [72]). For binaries with high mass ratio where a more massive binary also is more compact, this fraction 
approaches 30% for encounters with a periastron within ~ 20(fli + 02), essentially 3.5 times more efficient than equal 
mass binaries. For binary fractions of a few per cent, triple formation for a binary with a BH is more frequent than 
strong encounters with single stars [47]. It has been shown that for a seed BH-WD binary with an initial separation of 
a < SORq, the triple formation with subsequent Kozai cycle coupled with gravitational wave emission will result in an 
LMXB formation during several Gyr. 

CONCLUSIONS 

In studies of binary fractions, it is crucial to distinguish between the "observable" and the "true" binary fractions, they 
have different values and have different evolution. Soft binaries can not be ignored in simulations: they act as an energy 
sink and could be the key in achieving steady binary-burning phase. Triples play very important role in the formation 
of compact binaries, specifically. X-ray binaries with a NS or a BH accretor. 
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